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PESIOME 


Ареал представителей комплекса Darevskia (caucasica) (D. alpina, Р. с. caucasica, Р. с. vedenica, Р. daghes- 
tanica) охватывает Главный Кавказский хребет и ero отроги на южном и северном макросклонах B пре- 
делах высот от 500 м до 3200 м над уровнем моря. В работе проанализирован участок митохондриаль- 
ного гена цитохрома b (1096 п.н.) от 48 особей этих ящериц, собранных в 23 локалитетах. Обнаружены 
39 уникальных гаплотипов. Минимальные значения показателей генетического разнообразия выяв- 
лены у Р. с. vedenica, максимальные — у D. с. саисаяса. Все таксоны комплекса, за исключением ве- 
денской ящерицы, демонстрируют генетическую гетерогенность. На дендрограмме филогенетических 
отношений можно выделить четыре клады, соответствующие D. alpina, D. c. caucasica, D. c. vedenica, 
р. daghestanica. Р. с. caucasica формирует две обособленные гаплогруппы. Выборка веденской ящерицы 
сформировала самостоятельную кладу на дереве, близкую по положению к номинативному подвиду 
кавказской ящерицы. По митохондриальному маркеру комплекс Darevskia (caucasica) полифилетичен: 
нуклеотидные последовательности D. alpina сформировали единую кладу с Darevskia (saxicola), обра- 
зовав субкладу с последовательностями D. brauneri. Пример популяции D. daghestanica, обнаруженной 
в лесном биотопе без скальных выходов и генетически не отличающейся от особей с сопредельной тер- 
ритории, свидетельствует о способности скальных ящериц осваивать не свойственные им биотопы за 
короткий период времени. Генетическая дистанция между D. с. vedenica и D. с. саисаяса сопоставима 
с таковой для пары Р. с. caucasica и Р. daghestanica. Мы предполагаем наличие самостоятельной формы 
D. саисаяса ssp., обитающей в Южной Осетии и на сопредельных территориях. 


Ключевые слова: Кавказ, Darevskia alpina, Dareoskia caucasica, Darevskia caucasica vedenica, Darevskia 
daghestanica 
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ABSTRACT 


The range of superspecific complex of rock lizards Darevskia (caucasica) (D. alpina, D. c. caucasica, D. c. vedeni- 
ca, D. daghestanica) includes the Main Caucasian ridge and its spurs on the southern and northern macroslopes 
within heights from 500 m to 3200 m above sea level. A fragment of the mitochondrial cytochrome b gene (1096 bp) 
was analyzed in 48 specimens of these lizards. Specimens were collected in 23 localities. We have found 39 unique 
haplotypes. The minimum genetic diversity was found in D. c. vedenica, maximum in D. c. caucasica. All taxa of 
the complex, except for D. c. vedenica, demonstrate genetic heterogeneity. Four clades can be distinguished on the 
dendrogram of phylogenetic relationships. They correspond to D. alpina, D. c. caucasica, D. c. vedenica, D. dagh- 
estanica. D. c. caucasica forms two distinct haplogroups. On the phylogenetic tree, the sample of D. c. vedenica 
formed an independent clade close to D. c. caucasica. The Darevskia (caucasica) complex is polyphyletic accord- 
ing to the analysis of the mitochondrial marker. Nucleotide sequences of D. alpina formed a single clade with 
Darevskia (saxicola) complex representatives and combined with the sequences of D. brauneri into subclade. On 
the example of the population of D. daghestanica, found in a forest biotope without rocky outcrops and genetically 
not differing from individuals from the adjacent territory, the ability of rocky lizards to master biotopes unusual 
for them in a short period of time was demonstrated. This example demonstrates the ability of rock lizards to adapt 
to unusual biotopes in a short period of time. Genetic distance between D. c. vedenica and D. c. caucasica is com- 
parable to the distance between D. c. caucasica and D. daghestanica. We assume the existence of an independent 
form of D. caucasica ssp. inhabiting South Ossetia and adjacent territories. 


Key words: the Caucasus, Darevskia alpina, Darevskia caucasica, Darevskia caucasica vedenica, Darevskia 
daghestanica 





ВВЕДЕНИЕ 


Надвидовой комплекс Darevskia (caucasica), 
включающий альпийскую D. alpina (Darevsky, 
1967), кавказскую Р. с. caucasica (Méhely, 1909), 
веденскую D. с. vedenica (Darevsky et Roitberg, 
1999) и дагестанскую D. daghestanica (Darevsky, 
1967) ящериц (Puc. 1), можно рассматри- 
вать Kak сложную в систематическом плане и 


недостаточно изученную группу (Roytberg 1999; 
Лотиев и Доронин [Lotiev and Doronin] 2011; Ty- 
ниев и др. [Tuniyev et al.] 2017). Ареал этих так- 
сонов охватывает Главный Кавказский хребет 
иего отроги на южном и северном макросклонах 
в пределах высот от 500 м до 3200 м над уров- 
нем моря (Даревский [Darevsky] 1967; Доронин 
[Doronin] 2013a; данные авторов). Сложная оро- 
графия и гидросеть, гетерогенность климата 
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Рис. 1. Представители комплекса Darevskia (caucasica): А — D. alpina, самец, урочище Джилы-Су, Кабардино-Балкария 
(фото И.В. Доронина, 8.06.2019); В — D. daghestanica, самец, Аргунское ущелье, Чечня (фото К.Ю. Лотиева, 14.05.2018); 
С - D. с. caucasica, Ардонское ущелье, Северная Осетия-Алания (фото К.Д. Мильто, 14.08.2006); D - D. c. vedenica, самец, 
12 км южнее c. Рошни-Чу, Чечня (фото К.Ю. Лотиева, 09.05.2018); Е — D. с. caucasica, неполовозрелая особь, самки и самец, 
окрестности c. Бритат, Южная Осетия (фото Б.С. Туниева, 24.06.2014). 


Fig. 1. Representatives of the Darevskia (caucasica) complex: А — D. alpina, male, Dzhily-Su, Kabardino-Balkaria (photo by I.V. Doro- 
nin, 8.06.2019); B — D. daghestanica, male, Argun Gorge, Chechnya (photo by K.Yu. Lotiev, 14.05.2018); C — D. c. caucasica, Ardon 
Gorge, North Ossetia-Alania (photo by K.D. Milto, 14.08.2006); D — D. c. vedenica, male, 12 km south of Roshni-Chu village, Chechnya 
(photo by K.Yu. Lotiev 09.05.2018); E — D. caucasica, juvenile, females and male, vicinity of Britat village, South Ossetia (photo by 
B.S. Tuniyev, 24.06.2014). 
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Рис. 2. Пункты сбора экземпляров Darevskia (caucasica), использованных для молекулярно-генетического анализа: 
D. alpina (красные кружки), D. с. caucasica (зеленые кружки), D. с. vedenica (синие кружки) и D. daghestanica (желтые круж- 
ки). Нумерация соответствует таковой в Табл. 1. 


Fig. 2. Geographic distribution of the samples of Darevskia (caucasica) used in the molecular-phylogenetic analysis: D. alpina (red circ- 
les), D c. caucasica (green circles), D. c. vedenica (blue circles) and D. daghestanica (yellow circles). Numbers of localities are identical 


to those in the Table 1. 


и растительного покрова, создающие большую 
мозаичность ландшафтов этого горного ре- 
гиона, а также его насыщенная палеоистория 
предполагают сложную филогеографическую 
картину у обозначенного комплекса. На это 
указывает и вероятный реликтово-островной 
характер распространения скальных ящериц 
(Tuniyev et al. 2019). 

Остается открытым вопрос вхождения Р. alpi- 
па в данный комплекс. В многочисленных публи- 
кациях по филогении скальных ящериц (Fu et al. 
1995, 1997; Fu 1999; Murphy et al. 2000) была ис- 
пользована выборка альтийской ящериц с горы 
Аишхо (Краснодарский край, Россия), которая 
по результатам анализа мтДНК вошла в «saxi- 
cola group». Д.Г. Чобану с соавторами (Ciobanu 
et al. 2003), изучавшие сателлитную ДНК, также 
включили этот вид в кладу «зам со а». 

За последнее время вышла только одна пу- 
бликация, посвященная применению молеку- 


лярно-генетических методов для изучения фи- 
логенетических взаимоотношений в комплексе 
Darevskia (caucasica) (Доронин и др. [Doronin et 
al.] 2016). Работа показала эффективность ДНК- 
баркодинга (анализа на основе изменчивости 
фрагмента митохондриального гена субъедини- 
цы І цитохром оксидазы) для идентификации 
D. alpina, D. c. caucasica и D. daghestanica; была 
продемонстрирована значительная дифферен- 
циация внутри D. alpina по сравнению с даге- 
станской и кавказской ящерицами. Однаковвиду 
ограниченности изученного материала (24 экз. 
из 9 локалитетов) исследование не охватило всех 
представителей комплекса. Это касается послед- 
него по времени описания таксона — Р. с. vedenica, 
статус которого ставился под сомнение (Кузьмин 
и Семенов [Kuzmin and Semenov] 2006). Кроме 
того, новые находки веденской ящерицы в преде- 
лах Чечни (Лотиев [Lotiev] 2020), D. daghestanica 
на севере Азербайджана (Doronin et al. 2019; 
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Мустафаева и др. | Mustafaeva et al.] 2019), а также 
полученные сведения о морфологических отли- 
чиях D. caucasica с территории Южной Осетии от 
северокавказских представителей вида (Туниев 
и др. [Tuniyev et al.] 2017) ставят новые вопросы B 
изучении комплекса. Таким образом, назрела не- 
обходимость более обстоятельного изучения фи- 
логении и филогеографии ящериц этой группы. 
Полученные сведения могут быть исполь- 
зованы как при реконструкции путей и эта- 
пов фауногенеза региона, так и при планиро- 
вании природоохранной деятельности. Уже 
сейчас ряд таксонов комплекса занесены в ре- 
гиональные и национальные Красные книги 
(Доронин [Doronin] 2013; Туниев и Туниев 
[Tuniyev and Tuniyev] 2012, 2017; Лотиев и Ту- 
Hues [Lotiev and Tuniyev] 2017; Лотиев [Lotiev] 
2020). Молекулярно-генетические данные по- 
зволяют выделить наиболее уникальные попу- 
ляции и обосновать необходимость их охраны 


(González-del-Pliego et al. 2019). 


МАТЕРИАЛ И МЕТОДЫ 


В качестве образцов тканей ящериц исполь- 
зовали части аутотомированного хвоста или 
фаланги пальцев передних конечностей, фикси- 
рованные в 96% этаноле, от 47 экз., собранных 
в 22 локалитетах (Рис. 2). Полученные нами по- 
следовательности (MW583669—583715) были 
депонированы в GenBank NCBI (http://www. 
ncbi.nlm.nih.gov/gene/). Ббльшая часть Bayuep- 
ных экземпляров хранится в коллекции 300- 
логического института РАН (ZISP) и Инсти- 
тута зоологии Академии Наук Азербайджана 
(IZA) (Табл. 1). Кроме того, была использована 
последовательность D. саисаяса из GenBank — 
MH481131 (Tarkhnishvili et al. 2020). В ходе моле- 
кулярно-генетического анализа изучены после- 
довательности фрагмента полноразмерного гена 
цитохрома b мтДНК (cyt b), широко применяе- 
мого в таксономических исследованиях ящериц 
семейства Lacertidae, в том числе и рода Darevskia 
Arribas, 1999! (Доронин и др. [Doronin et al.] 2013; 
Freitas et al. 2016; Kurnaz et al. 2019 и др.). 





! Согласно решению Международной комиссии IIO зоологи- 


ческой номенклатуре датой опубликования описания рода 
Darevskia следует считать не 1997, a 1999 г. (International 
Commission on Zoological Nomenclature 2020). 


Выделение тотальной ДНК из образцов TKa- 
ней проводили с помощью коммерческого набо- 
pa DNeasy Blood & Tissue Kit (QIAGEN) соглас- 
HO протоколу производителя. Амплификацию 
участка гена cyt b проводили с использовани- 
ем T100TM Thermal Cycler (Bio-Rad) и прайме- 
ров 5’ - AAC CRC УСТ TGT MTT CAA CTA - 3’ 
(LgLu) и 5’ - GGY TTA CAA GAC CAG YGC 
ТТТ - 3’ (RtHr) (Доронин и др. [Doronin et al.] 
2013) по следующей схеме: предварительная де- 
натурация 94°С (5 мин), затем 34 цикла (94°С — 
30c,50'C – 60 c, 72°C – 60 с) и конечная элонгация 
при 72°С (5 мин). Секвенирование проводили на 
генетическом анализаторе ABI PRIZM 3500xL 
(Applied Biosystems) с применением Tex же прай- 
меров, что использовали при амплификации. 
Окончательная длина амплифицированного 
фрагмента ДНК составила 1096 п.н. Выделение 
ДНК и амплификацию фрагментов проводили 
B pecypcHom центре СПбГУ «Хромас». Для вы- 
полнения секвенирования амплифицирован- 
ных фрагментов ПЦР-продукты передавали 
B ресурсный центр СПбГУ «Развитие молеку- 
лярных и клеточных технологий». 

Выравнивание последовательностей произ- 
водили в программах Geneious Prime 2021.0.1 
(https://www.geneious.com) и AliView 1.6 (Larsson 
2014). Филогенетические деревья были pe- 
конструированы в процессе Байесова анализа 
(Bayesian analysis) в программе Mr. Bayes 3.1.2 
(Huelsenbeck and Ronquist 2001; Ronquist and 
Huelsenbeck 2003; Altekar et al. 2004) u mero- 
дом максимального правдоподобия (Maximum 
Likelihood, ML) в программном пакете MEGA X 
(Kumar et al. 2018). Статистическую надежность 
узлов МГ-деревьев оценивали путем бутстреп 
анализа (1500 псевдореплик для каждого дере- 
ва). Выбор модели эволюции для нуклеотидных 
последовательностей проводили в программе 
jmodeltest2 (Darriba et al. 2012) с использовани- 
em Байесовского информационного критерия 
(BIC) и информационного критерия Акаике 
(А1С). Была выбрана общая модель обратно- 
го хода времени (General Time Reversible, GTR) 
с параметрами I = 0.3057 и С = 0.6292. 

Для получения графических изображений 
деревьев применяли программу FigTree 1.4.4 
(http: //tree.bio.ed.ac.uk/software/figtree). Для уко- 
ренения были выбраны последовательности 
живородящей ящерицы Zootoca vivipara (Lich- 
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Таблица 2. Генетическое разнообразие в выборках Darevskia (caucasica). 


Table 2. Genetic diversity in the samples of Darevskia (caucasica). 
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(Taxon) n H h+SD n+SD S/n K Fu's Fs 
D. alpina 7 6 0.952 + 0.096 0.0235 + 0.0049 63/63 25.81 1.76 
D. daghestanica 21 17 0.976 0.023 0.01176+ 0.0010 48/48 12.88 -3.35 
D. c. caucasica 12 10 0.970 + 0.044 0.0325 + 0.0031 77/79 35.69 1.48 
D. c. vedenica 8 6 0.929 + 0.084 0.0013 + 0.0002 5/5 1.42 -3.58 


Примечание: n — объем выборки; Н — число гаплотипов; h — разнообразие гаплотипов; л — разнообразие нуклеотидов 
(на сайт); S — общее число полиморфных позиций; N — общее число мутаций; К — среднее количество нуклеотидных замен; 


Fu's Fs - значение Fs-recra. 


Note: n -sample size; Н — number of haplotypes; h — haplotype diversity; x — nucleotide diversity (per site); S — total number of poly- 


morphic positions; п — total number of mutations; К — average number of nucleotide substitutions; Fu's Fs — Fs test value. 


tentein, 1823) (NC026867) и тянь-шаньской 
ящурки  Eremias  stummeri Wettstein, 1940 
(NC029878), депонированные в GenBank NCBI. 
При проверке гипотезы о полифилетичности 
комплекса были использованы сиквенсы пред- 
ставителей Darevskia (saxicola): ящерицы Бра- 
унера D. brauneri (Méhely, 1909) (]Х0461614— 
0461620) и скальной ящерицы D. saxicola 
(Eversmann, 1834) (AF206180, MG704921), cum- 
патричных с таксонами комплекса Darevskia 
(caucasica). 

При анализе генетической изменчивости 
с использованием программы DnaSP у.5.10.01 
(Librado and Rozas 2009) вычисляли следующие 
параметры: общее количество полиморфных 
позиций (S), общее количество замен (п), ко- 
личество гаплотипов (Н), разнообразие гапло- 
типов (h), разнообразие нуклеотидов на сайт 
(п), среднее количество нуклеотидных замен 
(К), значение Еѕ-статистики Фу (Киз Fs). Io- 
строение сети гаплотипов осуществляли мето- 
дом TCS в программе PopART (Leigh and Bryant 
2015) Ha основе оригинальных данных авторов. 
Генетические дистанции (р-дистанции) вычис- 
ляли B программе MEGA X. 


РЕЗУЛЬТАТЫ 


Среди 48 исследованных последовательно- 
стей гена cyt b представителей комплекса об- 
наружены 39 гаплотипов (Табл. 2). Выявлены 
285 вариабельных позиций (26% от общей дли- 
ны фрагмента), из них — 259 парсимониинфор- 
мативных (23.6%). Соотношение транзиций 
к трансверсиям составило 5.21. Нуклеотидный 


состав участка Ha Г-цепи смещен в сторону 
уменьшения доли гуанина (A — 26.6, T - 31.9, С — 
29.8, G — 11.7%). 

Показатели генетического разнообразия 
представителей комплекса приведены в Табл. 2. 
Минимальные значения выявлены у Р. с. veden- 
ica, а максимальные - у D. с. саисаяса. Показа- 
тель нуклеотидного разнообразия (л) довольно 
низок у всех таксонов, что может быть объясне- 
но небольшим размером некоторых выборок. 
Выборки D. c. vedenica и D. daghestanica имеют 
слабо отрицательные значения Fs-recra Фу. 

Средняя величина некорректированных ге- 
нетических дистанций (p-distance) между все- 
ми изученными выборками составила 11.13 + 
2.53%, минимальная — между D. c. саисаяса 
и Р. daghestanica (8.31%), максимальная — меж- 
ду Р. с. vedenica и О. alpina (13.98%) (Табл. 3). 


Таблица 3. Средние генетические р-дистанции (%) (под ди- 
агональю) и их стандартные отклонения (над диагональю) 
между кладами комплекса Darevskia (caucasica) по данным 
анализа последовательностей фрагмента гена cyt b. 


Table 3. The average genetic p-distances (%) (below the diagonal) 
and their standard deviations (above diagonal) between the 
clades of the Darevskia (caucasica) complex based on the analysis 
of the sequences of cyt b gene fragment. 





Таксон 





(Тахоп) 5 B © р 
А (D. c. vedenica) 0.75 0.87 1.01 
B (D. c. caucasica) 8.48 0.71 0.88 
С (D.daghestanica) 10.10 8.31 0.97 
D (D. alpina) 13.98 12.17 13.71 
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TS2957 D c vedenica Chechnya Roshni Chu 
TS2958 D c vedenica Chechnya Roshni Chu 
^ TS2960 D c vedenica Chechnya Roshni Chu 

TS2959 D c vedenica Chechnya Kharachoy Ж A D p C ve d e n i ca 
1/0/100 TS2962 D c vedenica Chechnya Fortanga Riv 
TS2956 D c vedenica Chechnya Kharachoy Ж 
ZISP30997 D c vedenica Chechnya Kharachoy Ж 
TS2954 D c caucasica S Osetia Zgubir 
TS2950 D c caucasica S Osetia Verkhiy Ruk 
TS2965 D c caucasica S Osetia Verkhiy Ruk 
TS2953 D c caucasica N Osetia Kassar gorge 
ZISP26315 D c caucasica Kab Balk Chegem gorge Н 
ZISP26319_D_c_caucasica_Kab_Balk_Chegem_gorge B D. C. caucasica 
ZISP26317_D_c_caucasica_Kab_Balk_Chegem_gorge 
MH481131_D_c_caucasica_Georgia_Krestovyi_pass * 
ZISP26316 D c caucasica Kab Balk Chegem gorge 
ZISP26600 D c caucasica Kab Balk Chegem gorge 
TS2942 D daghestanica Dagestan V Gakvari 
TS2941 D daghestanica Dagestan V Gakvari 
TS2940 D daghestanica Dagestan V Gakvari 
TS2947 D daghestanica Dagestan V Gakvari 
~ IZA1023 D daghestanica Azerb Afurja 
1ZA1099 D daghestanica Azerb Laza 
IZA1101 D daghestanica Azerb Laza 
1ZA1094 D daghestanica Azerb Yeddi Gozel П 
TS2963_D_daghestanica_Chechnya_Argun_gorge C D Li dag h estan 1 са 
TS2935_D_daghestanica_Dagestan_Shaitli 
TS2936_D_daghestanica_Dagestan_Shaitli 
TS2948_D_daghestanica_Dagestan_Shaitli 
ZISP22218 5 D daghestanica Dagestan Kishcha 
TS2943 D daghestanica Dagestan Chonkatau Ж 
TS2944 D daghestanica Dagestan Chonkatau Ж 
TS2945 D daghestanica Dagestan Termenlik 
TS2946 D daghestanica Dagestan Termenlik 
ZISP21607 2 D alpina Adygea Tybga 
ZISP27945 D alpina Kar Cher Sofia waterfalls 
TS2927 D alpina Kab Balk Terskol * D D a 1, р i n a 
ZISP27943 D alpina Kab Balk Dzhily Su " 
ZISP26548 D alpina Kar Cher Zakan 
ZISP26540 D alpina KrasnodarR Azish Tau 
NC 026867 Zootoca vivipara 





NC 029878 Eremias stummeri 





0.07 


Рис. 3. Дендрограмма филогенетических отношений представителей комплекса Darevskia (caucasica) по результатам 
анализа фрагмента гена cyt b (1096 п.н.). Для узлов указаны байесовы апостериорные вероятности и бутстрэп-поддержки 
(BA/BS). Звездочкой отмечены нуклеотидные последовательности особей, собранных максимально близко к типовому 
местонахождению таксона. Буквенное обозначение клад соответствует таковой в Табл. 3. 


Fig. 3. Dendrogram of phylogenetic relationships of Dareoskia (caucasica) based on cytochrome b gene sequencing. Bayesian posterior 
probabilities and bootstrap values (BA/BS) are given. The asterisk marks the sequences of the topotypes, or specimens, collected as 
close as possible to the type locality of the taxon. The letter designation of the clades corresponds to that in Table 3. 


По результатам реконструкции филогенети- точного Кавказа; кладу В составила выборка 
ческих взаимоотношений среди скальных яще- D.c. caucasica c Центрального Кавказа, в которой 
PHIL комплекса можно выделить четыре клады можно выделить две субклады (р-дистанция = 
(Рис. 3). Кладу А образует D. c. vedenica Boc- 5.36%); кладу C — восточнокавказская D. daghe- 
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1 sample 
e D. alpina 
© D. daghestanica 
Өр. с. саисаѕіса 
e» c. vedenica 











Рис. 4. Медианная сеть гаплотипов цитохрома b представителей комплекса Darevskia (caucasica) (1096 п.н.). Каждый цвет- 
ной кружок — уникальный гаплотип, его размер пропорционален встречаемости в выборке; цвет кружка соответствует 
видовой принадлежности. Соединительные линии — вероятные эволюционные связи, засечки (или цифры в скобках) 
Ha линиях — нуклеотидные замены. Черные кружки в узлах линий — предсказанные гаплотипы. Условные обозначения: 
А - D. c. vedenica, В — D. с. caucasica, C – D. daghestanica, D — D. alpina. Буквенное обозначение групп гаплотипов соответству- 
ет таковому на Рис. 3. 


Fig. 4. The median network of cytochrome b haplotypes in representatives of the Darevskia (caucasica) complex (1096 bp). On the 
network: each colored circle is a unique haplotype; its size is proportional to its occurrence in the material; the color of the circle corre- 
sponds to the species. The connecting lines are probable evolutionary links; the nicks (or numbers in brackets) on lines are the nucleotide 
substitutions. Black circles at the nodes of the lines are the predicted haplotypes. Legend: A — D. c. vedenica, B — D. с. caucasica, 


C - D. daghestanica, D — D. alpina. The letter designation of the haplotype groups corresponds to that in Fig. 3. 


stanica; и кладу D — западнокавказская и цен- 
трально-западнокавказская D. alpina, образую- 
щие две субклады (р-дистанция = 3.58%). 

На полученной дендрограмме (Рис. 3) образ- 
цы D. с. оедетса достоверно сформировали ca- 
мостоятельную кладу, близкую по положению 
к номинативному подвиду кавказской яще- 
рицы, при этом все таксоны комплекса, за ис- 
ключением веденской ящерицы, демонстриру- 
ют генетическую гетерогенность. Аналогичное 
прослеживается и на медианной сети гаплоти- 
пов (Рис. 4). D. с. саисаяса представлена двумя 
гаплогруппами, обособленными друг от друга 


на 52 нуклеотидные замены. Это соотносится 
с подразделением этого таксона на две субклады 
на филогенетическом дереве (Рис. 3). Гаплоти- 
пы, соответствующие D. с. vedenica (обозначены 
буквой А), отстоят от гаплотипов Р. с. caucasica 
(B) на 85 замен. Генетическая р-дистанция MEX- 
ду гаплогруппами В2 и А составляет 8.30%, 
а между B1 и A - 8.73%. D. daghestanica пред- 
ставлена тремя гаплогруппами, незначительно 
обособленными друг от друга (6—10 нуклеотид- 
ных замен), что хорошо соотносится с результа- 
тами, полученными Ha филогенетическом дере- 
ве. Последовательности альпийской ящерицы 
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Clade А D. c. vedenica 

Clade B D. c. caucasica 

Clade C D. daghestanica 
ZISP26548 D alpina Kar Cher Zakan 

ZISP26540 D alpina KrasnodarR Azish Tau 

ZISP21607 2 D alpina Adygea Tybga H 
ZISP27945 D alpina Kar Cher Sofia waterfalls D. alpina 
TS2927 D alpina Kab Balk Terskol 

ZISP27943 D alpina Kab Balk Dzhily Su 
JX041618 1 Darevskia brauneri 

JX041620 1 Darevskia brauneri 

JX041617 1 Darevskia brauneri 

JX041619 1 Darevskia brauneri 

JX041614 1 Darevskia brauneri 

JX041615 1 Darevskia brauneri 

JX041616 1 Darevskia brauneri 

AF206180 1 Darevskia saxicola 

MG704921 1 Darevskia saxicola 

NC 026867 Zootoca vivipara 


NC. 029878 Eremias stummeri 


Puc. 5. Дендрограмма филогенетических отношений представителей комплекса Darevskia (caucasica) и Darevskia (saxicola) 
по результатам анализа фрагмента гена cyt b мтДНК. Указаны байесовы апостериорные вероятности и бутстрэп-поддерж- 


ки (BA/BS). 


Fig. 5. Dendrogram phylogenetic relationships of Darevskia (caucasica) and Darevskia (saxicola) based оп cyt b gene sequencing. 
Bayesian posterior probabilities and bootstrap values (BA/BS) are given. 


подразделяются Ha две гаплогруппы (D1 и D2), 
отделенные друг от друга 23 нуклеотидными 
заменами. Гаплогруппы Р. alpina также совпа- 
дают с субкладами на филогенетическом дереве 
(Рис. 3). 

С целью проверки гипотезы о монофилетич- 
ности комплекса была построена дендрограмма 
филогенетических отношений с включением 
представителей Dareoskia (saxicola) (Рис. 5). Как 
ивпредшествующих публикациях по этому во- 
просу, изученные нуклеотидные последователь- 
ности D. alpina сформировали единую кладу 
с Darevskia (saxicola), образовав субкладу с по- 
следовательностями D. brauneri (р-дистанция = 
5.34%, что значительно меньше, чем между так- 
сонами Darevskia (caucasica)). 


ОБСУЖДЕНИЕ 


Наша диагностика изученных экземпля- 
ров на основе внешней морфологии (использо- 
ваны признаки из: Даревский [Darevsky] 1967; 


Darevsky апа Roitberg 1999) полностью совпала 
с их распределением по кладам на полученной 
дендрограмме филогенетических отношений 
и IIO гаплогруппам на сети гаплотипов. 

Топология дендрограммы в целом согласует- 
ся со схемой филогенетических связей в работе 
Доронина с соавторами ([Doronin et al.] 2016): 
подтверждается обособление основных клад 
и видовой уровень дифференциации D. cauca- 
sica и D. daghestanica. 

Какбыло показано выше, по результатам ана- 
лиза cyt b D. арта входит в комплекс Darevskia 
(saxicola), а не Darevskia (caucasica), что полно- 
стью соответствует выводам предшествующих 
исследователей. Этого нельзя сказать на основе 
данных по внешней морфологии: альпийская 
ящерица была описана как подвид кавказской 
и долгое время рассматривалась в этом стату- 
се (Darewskij 1984). Полученная картина может 
быть объяснена интрогрессией митохондри- 
альной ДНК вследствие гибридизации с сим- 
патричными видами, возникающей при слабой 
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Рис. 6. Распространение гаплогрупп cyt b мтДНК представителей комплекса Darevskia (caucasica): D. арта (красный цвет), 
Dc. caucasica (зеленый цвет), D. c. vedenica (синий цвет) и D. daghestanica (желтый цвет). Номера локалитетов соответствуют 
таковым в Табл. 1. 


Fig. 6. Distribution of cyt b mtDNA haplogroups of Darevskia (caucasica): D. alpina (red), D c. caucasica (green), D. c. vedenica (blue) 
and D. daghestanica (yellow). Numbers of localities are identical to those in the Table 1. 


репродуктивной изоляции у скальных ящериц, 
на что впервые указал Даревский ([Darevsky] 
1967). Показательно, что гадюка Динника Vi- 
pera dinniki (Nikolsky, 1913), схожая с D. alpina 
очертаниями ареала и типом распростране- 
ния в пределах горно-лугового пояса Большого 
Кавказа (Tuniyev et al. 2019), входит в сложный 
гибридный комплекс горных гадюк (Zinenko 
et al. 2015). Также нельзя исключать и феноти- 
пическую конвергенцию, известную для раз- 
ных групп скальных ящериц (Tarkhnishvili et al. 
2020). До проведения исследований с использо- 
ванием ядерных маркеров мы продолжаем рас- 
сматривать альпийскую ящерицу как предста- 
вителя комплекса Darevskia (caucasica). 
Несмотря на значительное географическое 
расстояние между изученными популяция- 
ми D. с. vedenica и О. daghestanica, достигающее 
60 км для первой и 325 км для второй, и предпо- 
лагаемый дизъюнктивный характер ареала ве- 
денской ящерицы (Лотиев и Доронин [Lotiev and 


Doronin| 2011), их выборки продемонстрировали 
низкие показатели генетического разнообразия. 
Это говорит о вероятном прохождении популя- 
ций через «бутылочное горлышко», которое про- 
исходит при резком сокращении численности 
и ареала. Согласно данным по генетическому 
разнообразию (Табл. 2) О. alpina и Р. с. caucasica, 
напротив, в недавнем эволюционном прошлом 
претерпели рост численности и расширение 
ареалов. В то же время для дагестанской ящери- 
цы мы наблюдаем наличие структурированно- 
сти – обособление гаплогрупп, локализованных 
на северных отрогах восточной гряды Главного 
Кавказского хребта в юго-восточной части Даге- 
стана (C1), южном склоне Главного Кавказского 
хребта на севере Азербайджана (С2) и на терри- 
тории от Аргунского ущелья Чечни на западе до 
предгорий Дагестана на востоке (СЗ) (Рис. 6). 
При этом их дистанцированность друг от дру- 
га значительно уступает таковой у гаплогрупп 
D. арта и D. с. caucasica. 
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Рис. 7. Биотоп D. daghestanica в окрестностях с. Манасаул, Дагестан (фото Л.Ф. Мазанаевой, 15.04.2012). 
Fig. 7. Biotope of D. daghestanica in the vicinity of the Manasaul village, Dagestan (photo by L.F. Mazanaeva, 15.04.2012). 


Достаточно неожиданным стало географи- 
ческое распределение выборки альпийской 
ящерицы: особи с хребта Азиш-Тау и окресно- 
стей пос. Закан — удаленных друг от друга лока- 
литетов, образовали гаплогруппу D2, тогда как 
особь с расположенной относительной близко 
от них горы Тыбга вошла в гаплогруппу D1 вме- 
сте с ящерицами из района Архыза и Приэль- 
брусья. Схожую картину расхождения выбор- 
ки, отражающую сложную историю расселения 
скальных ящериц, мы наблюдаем у D. с. caucasi- 
са (Рис. 6): особи из близко расположенных ло- 
калитетов в Грузии и Южной Осетии подразде- 
лились на гаплогруппы B1 и B2. 

В данном исследовании были задейство- 
ваны экземпляры D. daghestanica, собранные 
в окрестностях с. Манасаул Дагестана и возле 
туристической базы «Терменлик». Эта изолиро- 
ванная популяция представляет значительный 
интерес, так как вид обитает здесь в лесном био- 


топе без скальных выходов (Рис. 7), где он впер- 
вые был найден в мае 2012 г. В июле 2020 г. мы 
обнаружили их повторно в лесу на территории 
самой базы, a также в развалинах старых зданий 
на железобетонных плитах. Локалитет находит- 
ся в лесной зоне предгорного Дагестана на се- 
верном макросклоне Гимринского хребта, отно- 
сящегося к передовым хребтам Внешнегорного 
района Дагестана. Ранее лесные массивы в этом 
районе покрывали значительные площади, но 
в настоящее время их ббльшая часть вырубле- 
на (Добрынин [Dobrynin] 1925). Ящерицы были 
отмечены на стволах и под корой дуба грузин- 
ского (Quercus iberica) и черешчатого (О. robus), 
а также граба кавказского (Carpinus caucasica), 
поросшего мхом аномодоном утонченным (Ano- 
modon attenuates), и под валежником. Ha yuer- 
HOM маршруте фиксировали до 4 экз./200 м. 

Ряд видов скальных ящериц населяет и лес- 
ные биотопы (Даревский [Darevsky] 1967; 
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Рис. 8. Биотоп D. daghestanica в Аргунском ущелье, Чечня (фото К.Ю. Лотиева, 12.07.2011). 
Fig. 8. Biotope of D. daghestanica in the Argun gorge, Chechnya (photo Ьу K.Yu. Lotiev, 12.07.2011). 


Орлова [Orlova] 1975), Ho для данного вида по- 
добные местообитания ранее не были отмечены. 
И.С. Даревский ([Darevsky] 1967) ornec D. daghe- 
stanica к группе «сухолюбивых» видов, Hace- 
ляющих сухие и умеренно сухие скалы (обыч- 
но сложенные из осадочных горных пород), их 
подножия и склоны с кустарниковой и травя- 
нистой растительностью, выходы твердых (не- 
редко вулканического происхождения) пород и, 
реже, глинистые обрывы в горно-степной и гор- 
но-луговой зонах (Рис. 8). Вероятно, дагестан- 
ская ящерица могла проникнуть в лесной пояс 
Гимринского хребта с Левашинского плато по 
склонам ущелий рек Пулеувла и Губденозень 
(хребет Чонкатау). На это указывает близость 
по молекулярно-генетическим данным особей 
из района Терменлик и хребта Чонкатау. Этот 
пример говорит о способности скальных яще- 
риц осваивать не свойственные им биотопы за 
короткий период времени. 


Обращает на себя внимание высокий уро- 
вень генетической дистанции между веден- 
ской и номинативной кавказской ящерицами 
(8.48%), сопоставимый с таковым для пары но- 
минативной кавказской и дагестанской ящериц 
(8.31%). Эти показатели превышают межви- 
довые дистанции для многих представителей 
рода, к примеру в комплексе Darevskia (saxicola) 
(Доронин и др. [Doronin et al.] 2013) и Darevskia 
(parvula) (Kurnaz et al. 2019). Дальнейшие иссле- 
дования филогении комплекса, B TOM числес HC- 
пользованием ядерных маркеров, позволят сде- 
лать вывод об их таксономическом положении. 
Нельзя исключать придание D. с. vedenica видо- 
вого статуса. 

В дальнейшем особое внимание следует 
уделить вопросу взаимоотношения двух под- 
видов кавказской ящерицы в предполагаемой 
зоне симпатрии, которая в настоящее время 
неизвестна. Возможность контакта двух форм 
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подтверждается недавним обнаружением вы- 
сокогорной субальпийской популяции веден- 
ской ящерицы на северном склоне горы Дайхох 
в Чечне (Лотиев [Lotiev] 2020). Напомним, что 
именно изучение симпатрических популяций 
р. саисаяса и О. daghestanica с территории Чеч- 
ни позволило говорить об их видовом статусе 
(Roitberg, 1994, 1999; Roytberg Lotiev, 19927). 

Аналогична ситуация и с югоосетинскими 
популяциями D. caucasica: ранее уже было пока- 
зано ихзначительноесвоеобразие,проявляюще- 
еся не только в особенностях окраски и рисунка 
(в частности, доминирование зеленоокрашен- 
ных самцов и наличие зеленоватых самок), но 
ив характере фолидоза (Рис. 1Е) (Туниев идр. 
[Tuniyev et al| 2017). Согласно полученным 
нами данным их представители сформировали 
самостоятельную субкладу (гаплогруппа В1 на 
Рис. 4) совместно с особью из Кассарского уще- 
лья в Северной Осетии, географически близко- 
го ктерритории Южной Осетии; она значитель- 
но отличается от северокавказской субклады 
(гаплогруппа B2 на Рис. 4) —р-дистанция между 
ними равняется 5.36%. Вероятно, в данном слу- 
чае речь может идти о новом подвиде. 

О выделении новых подвидов в <Lacerta 
caucasica group» писал Даревский (Darevsky 
1993), включивший в эту группу L. с. caucasica, 
«L. c. subsp. 1» и «L. c. subsp. 2». Под одним из этих 
подвидов, вероятнее всего, подразумевалась 
описанная через шесть лет веденская ящерица, 
а под другим — закавказские популяции кавказ- 
ской ящерицы. 

Отметим, что особи D. caucasica, схожие по 
окраске и элементам фолидоза с югоосетински- 
ми, обнаружены в окрестностях с. Камилух на 
р. Джурмут (приток р. Андийское Койсу, Даге- 
стан) (Туниев и др. | Tuniyev et al.| 2017), а нами 
были выявлены на сопредельной территории 
Грузии — вдолинер. Арагви (Мцхета-Мтианети). 

Не меньший интерес представляет изучение 
сим- и парапатричных популяций D. caucasica 
и D. daghestanica высокогорий Дагестана, где 
морфологическая разобщенность этих видов 
выражена довольно слабо, что может свидетель- 
ствовать о гибридизации между ними (Roitberg 


1999). 





? Вэтих публикациях речь идет о не описанной на тот момент 


D. c. vedenica. 


Сопоставление клад и гаплогрупп (Puc. 3) 

с типовыми местонахождениями таксонов ком- 

плекса (Доронин [Doronin] 2014) показало сле- 

дующее: 

— D. alpina описана из окрестностей с. Терскол 
Кабардино-Балкарии (43°15'N 42°30"E, по 
месту сбора голотипа); последовательности 
топотипов (ZISP TS 2927, TS 2928) вошли 
в субкладу, соотносящуюся с гаплогруппой 
D1, вместе с ящерицами из Адыгеи и Карача- 
ево-Черкесии; 

— D. caucasica описана с горы Казбек в Грузии 
(= 42°40'N 44°35"E, по месту сбора лектоти- 
па); последовательность от наиболее близко 
собранного к этому локалитету экземпля- 
ра с южного склона Крестового перевала 
(MH481131), ведущего из долины р. Терек 
в долину р. Арагви, вошел в субкладу, соот- 
носящуюся c гаплогруппой B2, вместе c яще- 
рицами из Кабардино-Балкарии; 

— D. с. vedenica описана по экземпляру (голо- 
типу), пойманному у дороги Ведено-Харачой 
(= Хорачой) в 1 км южнее с. Ведено Чечни 
(42°57'М 46'07'E) последовательности от 
наиболее близко собранных к этому лока- 
литету экземпляров из с. Харачой (ZISP 
TS 2956, TS 2959, ZISP 30997), как и все 
имеющиеся в нашем распоряжении после- 
довательности таксона, вошли в кладу А, не 
делящуюся на субклады; 

— D. daghestanica описана из c. Леваши (= Лава- 
ши) Дагестана (42°25'М 47°19'Е, по месту сбо- 
ра лектотипа); последовательности от наи- 
более близко собранных к этому локалитету 
экземпляров из Левашинского р-на (ZISP 
TS 2943, TS 2944) вошли в гаплогруппу C3 
вместе с ящерицами из Дагестана и Чечни. 
Таким образом предполагаемое описание но- 

вого подвида D. саисаяса ssp. с территории За- 

кавказья не нарушит стабильности номенкла- 
туры данного комплекса. 
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